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Recent developments on the rotational instabilities of relativistic stars are reviewed. 
The article provides an account of the theory of stellar instabilities with emphasis 
on the rotational ones. Special attention is being paid to the study of these insta- 
bilities in the general relativistic regime. Issues such as the existence relativistic 
r-modes, the existence of a continuous spectrum and the CFS instability of the 
w-modes are discussed in the second half of the article. 

1 Introduction 

The oscillations and instabilities of relativistic stars gained a lot of interest 
in the last decades because of the possible detection of their associated grav- 
itational waves. It is not impossible that every compact star in a specific 
period of its life will undergo an oscillatory phase, in which it becomes unsta- 
ble. Only the dynamical instabilities were thought to be a relevant source of 
gravitational waves, whereas instabilities due to dissipation mechanisms were 
believed to be only of academic interest. This belief has been dramatically 
altered during the last five years after it was discovepnd that for a specific 
class of rotational perturbatkms, the so called r-modeil^H, the instability due 
to gravitational radiationEHj'BH has the potential of being a prime source for 
gravitational waves. It was subsequently shown that this instability has many 
interesting astrophysical implications, which attracted the attention of both 
relativists and astrophysicists. 

In the first section of this review we shortly describe the various stellar 
instabilities and how they operate, for a more detailed review one can refer to 
a recent article by StergioulasQ. The second section is devoted to the r-mode 
instability with the focus on its application to gravitational-wave research arid 
astrophysics, for more details see a recent review by Andersson and Kokkotasu. 
The last section is devoted to an analysis of new features brought about by the 
relativistic treatment of the problem. Questions on the existence relativistic 
r-modes and a continuous spectrum are discussed in detail. A new result that 
is presented is the existence of the w-mode instability in ultra-compact stars, 
i.e. the instability due to the existence of an ergosphere of the spacetime or 
w-modesO. 
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2 Stability of Relativistic Stars 



2.1 Non-rotating stars 

There are many different types of instabilities whicfi can operate in a relativis- 
tic star. The most familiar instability is due to the existence of a maximum 
mass beyond which the star collapses. In other words it becomes unstable 
with respect to radial perturbations. The points of instability can be easily 
found by constructing a sequence of stellar models for a particular equation 
of state (EOS), parametrized for example by the central density Ej,- One can 
thus obtain the Mass- Radius diagram, where the extrema signal the transition 
from stability to instability and vice-versa. Within a perturbative approach, 
the stability analysis with respect toj^adial pulsations is reduced to a second- 
order Sturm-Liouville type equationllSlO with eigenvalue . For cr^ > the 
eigenfrequency a is real and the mode is stable while for cr^ < the eigen- 
frequency a becomes imaginary and the mode is unstable. The case — 
marks the onset of the instability. For the special case of uniform density 
stars, general relativity predicts that the star becomes unstable with respect 
to collapse if M/i? > 8/9. 

For non-radial perturbations of non-rotating stars a standard criterion for 
mode stability is the sign of the so called Schwarschild discriminant 



where r is the radial coordinate, p the pressure of the fluid and e the total 
energy density. The condition > everywhere inside the star is necessary 
for stability. If As < somewhem. inside the star, the g-modes become un- 
stable with respect to convectionEHl. The stability of non-radially oscillating 
stars can also be studied via perturbation theory. One cpi derive a system 
of perturbation equations for the fluid and the spacetimalj. For relativistic 
stars the frequency of oscillation modes is no longer real valued because the 
non-radial oscillations generate gravitational radiation which dissipates en- 
ergy from the star. Instead the frequency is complex [a = -\- isigmaj) and 
the imaginary part of the frequency describes the damping or the growing 
time of a mode. Unstable modes are characterized by t < 0. For non-rotating 
relativistic stars with As > no unstable modes have been found, although a 
rigorous proof of stability is still lacking. 
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(1) 



2.2 Rotating Stars 



The size of a rotating configuration is limited by the condition that the cen- 
trifugal acceleration should not exceed the centripetal acceleration of gravity. 
This leads to the following upper limit for the angulau^elocity in terms of 
the mean density e, established already a century agoEJ: 

< 27rGe. (2) 

The limiting rotation frequency, il.K, is called the Kepler limit. Although 
derived in Newtonian theory using simplifying assupuptions, this relation is 
approximately correct for relativistic systems as welEJ'ta. 

For uniformly rotating stars, stability against axisymmetric perturbations 
can be determined, via turning points along a constant angular momentum se- 
quence til. An initially uniformly rotating relativistic star under the influence 
of such an axisymmetric instability becomes differentially rotating. Since the 
viscosity is the main mechanism for the redistribution of angular momentum 
within the fluid, the timescale for the development of the instability is set by 
viscous dissipation. 

Rotational instabilities of stars arise from non-axisymmetric perturbations 
of the form e™"^, where 4> is the azimuthal angle. The bar mode (m = 2) is 
known to be the fastest growing unstable mode. There are two types of 
rotational instabilities: the dynamical and the secular ones. 

• The dynamical instability is driven by hydrodynamical and gravita- 
tional forces and develops on an extremely short timescale, the dynamical 
timescale. For neutron stars it is of the order of milliseconds. Of special 
interest is the so called bar mode instability {I = m = 2) which deforms the 
star into in a bar shaped structure within approximately one rotational 
period. The bar mode instability^ Ls an extremely efficient mechanism of 
producing gravitational wavest§llaE3. 

• The secular instabilities are driven by dissipative processes such 
as viscosity or gravitational radiation. They develop on much longer 
timescales, which are of the order of the dissipation time (for example 
the time required for viscosity to redistribute the angular momentum), 
i.e. tenths of seconds or even more. 

The viscosity driven instability plays hardly any role in neutron stars 
because it occurs only in|S±j(iLlar models with much stiffer EOS than those 
considered to be realistididS. 

ChandrasekharSB first proved that gravitational radiation can drive 
Maclaurin spheroids, that are rigidly rotating uniform density spheroids, 



florence: submitted to World Scientific on February 7, 2008 3 



unstable. This instability causes the Maclaurin spheroids to evolve into-a 
stationary but non-axisymmetric configuration, the Dedekind ellipsoidBj. 
This inijCpisting result was put on a rigorous footing by Friedman and 
SchutztrljtJ, who also proved that this instability is generic and that all 

rotating perfect fluid stars are unstable with respect to gravitational radi- 
ation. 



2.3 The CPS instability 

The mechanism far the gravitational-wave instability can be understood in 
the following wayOi3. Consider first a non-rotating star. Then the mode- 
problem leads to eigenvalues for cr^, i.e. one obtains equal values ±|cr| for 
the forward and backward propagating modes (corresponding to m — ±|m|). 
These two mode branches are affected by rotation in different ways. The 
retrograde mode will be dragged forward by the stellar rotation. If the star 
spins sufficiently fast, this mode will appear moving forward in the inertial 
frame, but still backward in the rotating frame. Thus, an observer at infinity 
sees gravitational waves with positive angular momentum emitted by this 
retrograde mode, but since the perturbed fluid rotates slower than it would 
in absence of the perturbation, the angular momentum of the mode itself is 
negative. The emission of gravitational waves consequently makes the angular 
momentum of the mode increasingly negative thus leading to an instability. 

From the above, one can easily conclude that a mode will be unstable if 
it is retrograde in the rotating frame (an observer rotating with the star and 
observing a frequency a for the mode) and prograde for a distant observer 
measuring a mode frequency a — mf2, i.e. if 

a{a - mn) < 0. (3) 

For the high frequency (/ and p) modes this is possible only for extremely 
high values of or for quite large m. In general, for every mode there will 
always be a specific value of m for which the mode will become unstable. 

This class of frame- dragging instabilities is usually referred to as 
Chandrasekhar-Friedman-Schutz (CFS) instabilities. It is easy to see that 
the CFS mechanism is not unique to gravitational radiation. Any radiative 
mechanism will have the same effect. 

A new result described in Section 4.3 is that the CFS mechanism is not 
only active for fluid modes but also for the spacetime or the so called w-modesEI. 
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Figure 1. A schematic illustration of the conditions under which the CFS instability is 
operating. A perturbed star can be viewed as a superposition of a uniformly rotating 
background and a non-axisymmetric perturbation. A mode is unstable if it is retrograde 
according to an observer rotating with the fluid (left), but appears prograde for an observer 
at infinity (right). 



2.4 Dissipative Effects 

Dissipation is important for the stability/instability of rotating stars. In gen- 
eral dissipation due to the fluid viscosity tends to suppress the instability 
and thus stabilizes stars that do not rotate fast enough. The typical fluid 
viscosities that have been used for the study of unstable / and r-modes are 
the shear and bulk viscosityE^. The effects of the dissipation on the evolution 
of a perturbation are most easily studied by considering the energy of the 
perturbation as measured in the rotating frame 



where 6u and 5$ are the perturbations of the fluid velocity and gravita- 
tional potential. In the absence of dissipation the energy of the mode is 
conserved. For a mode which is subject to dissipation, the eigenfunction is 
~ exp[z(fT -|- i/r)t], where a is the real part of the mode frequency and 1/t its 
imaginary part. The problem of estimating the stability of a mode is reduced 
to determining the sign of the imaginary part of the frequency. The imaginary 
part of the frequency is the ratio between the energy of the mode and the rate 
at which a specific dissipation mechanism dissipates energy from the mode, 
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(4) 



i.e. 



- = (5) 
T 2E dt ^ ' 

One now has to compare the damping/growth times due to the various mech- 
anisms for energy dissipation in a given neutron star model 

1111 

- - — + + . (6) 

T '^gw '^shcar 'H^ulk 

Rough estimates of the the various viscous mechanisms are given below. 

• Shear Viscosity is the main dissipation mechanism for low temperatures 
(T < l(fK) and it arises from the momentum transport due to neutron- 
neutron scattering. The dissipative energy loss in the mode energy is 
given by 

— = 2 y ^ Sa'^P Sa:^ dV, (7) 

where aaf3 is the shear tensor and the thermodynamic function rj repre- 
seiits-.the shear viscosity of the fluid. A typical value for neutron stars 



isl 



sE§E2l 



r, = 2x IQi^e^/Vg-^ g cm"! s'^. (8) 

Bulk Viscosity dominates at high temperatures (T > lO^K) when mat- 
ter becomes transparent to neutrinos. As the oscillating fluid undergoes 
compression and expansion, the weak interaction requires a relatively long 
time to re-establish equilibrium. This creates a phase lag between de 
sity and pressure perturbations which results to a large bulk viscosity! 
Typically the mode energy lost through bulk viscosity is radiated away 
by neutrinos. This energy loss can be estimated by: 

^ = 1 (6aSa* dV, (9) 

where Sa is the expansion associated with the mode and C the coefficient 
of the bulk viscosity, which for mode frequencies (in the rotating frame) 
has the value 

C = 6 X 1025(a + mn)-^^ el, g cm"! s-\ (10) 



6 
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• Gravitational Radiation dissipates energy from an oscillating star at a 
rate estimated via the post-Newtonian multipole-formulaES as 

^ = ^a{a + mn)Y,Ni<J^' (|<5A™|' + , (11) 



dt 

1=2 



where 



5 Jim = - ( 7^ ) I r' {eSu + u6e) ■ Yi^:dV, (14) 



_ 47rG (^ + l)(; + 2) 
' c2'+U(Z- 1)[(2? + 1)!!]2- ^ 

The radiation emitted is the sum of the mass 5 Dim and current 5 Jim 
multipole contributions given by 

5Dun = j 5tr' YCmdV, (13) 

1/2 

c\l + l 

where the -54^ are the spherical harmonics and Y^f^* the magnetic type 
multipolesEEIo. 

By estimating the energy of the mode using equation (^) and the energy 
dissipation via (0) , (^ and ( [Tl] ) one gets all the components needed in equation 
(^. For a given EOS and a given stellar model the sign, which signals if 
the mode is growing or not, of equation (^ is a function of the rotational 
period and the temperature of the star. Figure 2 shows the instability window, 
obtained in this way, of the I = m = 2 r-mode. 

3 r-modes 

The oscillation patterns of a star are classified according to the dominant 
restoring force that tries to push a displaced fluid element back into its equi- 
librium position. The two main families of modes, which exist for any type 
of stars (rotating or not), are those whose restoring forces are pressure {p- 
modes) and buoyancy (g-modes). Rotation not only shifts the spectra of 
these modes but additionally gives rise to a new type of restoring force, the 
Coriolis force, with an associated new family of rotational modes. In general 
these are called inertial modes and are known in geophysics as Rossby waves. 
Of special interest to us is the quadrupole r-mode with I — 2, m — 2. 

Inertial modes are primarily velocity perturbations. In the slow rotation 
limit, the velocity field of the r-mode is given by 

Ju^ai?f}(^)'l^fe-*, (15) 
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where a is the mode amphtude and R the radius of the star. In a first 
approximation the fluid elements have no radial displacement and move on 
paths which are approximately ellipses with ^-dependent eccentricities 

-\- + , ^ , ^ a'R^ (16) 
sin 9 sin 6 cos^ 6 

see for example Fig. 1 ini. The associated density perturbations are small 
for low rotation rates Sp ^ 0{^'^) and at the Newtonian level, the mode 
frequency in the rotating frame of reference is 

= WTTy 

Using the criterion for CFS stability given by equation we get 

a{o mil) = < 0, (18) 

which implies that the r-mode is unstable for any rotation rate of the star. 
This was first pointed out 4 years agoBH and immediately created considerable 
excitement because of the astrophysical implications this instability would 
have for the spin evolution of newly born neutron stars and for the search for 
gravitational waves. 

The main gravitational wave output from r-modes is not due to distur- 
bances in the star's fluid density (mass multipoles, SDim ~ ) but originates 
from the time dependent mass currents (current multipoles, SJim ^ This 
is the gravitational analogue of magnetic multipole radiation and is a rare case 
among the expected astrophysical sources of gravitational radiation. 

3.1 Instability window 

One can estimate the growth time of the instability using equation (||) together 
with the viscous damping times. In a first approximation one obtains 

Tgw ^ 47 M^l R^^^ P^^^ sec, (19) 

- 7 X 10^ Af-f * R^^'^ T2 sec, (20) 

Tbv 3 X 10" Mi.4 R^Q Pis Tg^ sec. (21) 

From the above values one can easily show that for high temperatures 
(T > 10^"K) the bulk viscosity suppresses any oscillations, while the shear 
viscosity dominates for temperatures below lO^K. Still, as shown in Fig. 2, 
for temperatures betwpejj_|_0'' and 10^*^, there is a broad window in which the 
instability can operat 
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Figure 2. The instability window for the I = m = 2 r-mode of a typical neutron star 
{R = 10 km and M = IAMq and Kepler period Pk ~ 0.8 ms). 



In particular these results imply that for temperatures of the order of 
10^ K a small rotation rate of the order of 5% of the Kepler frequency is 
already sufficient to drive the star unstable. 

The growth time due to the mass multipole radiation of the r-modes 
is about four orders of magnitude longer and the corresponding instability 
window is considerably smallem. The picture is similar for r-modes with 
higher I. 

In order to obtain a more realistic description of neutron stars, one should 
also take into account the effects of a solid crust. In old neuirnji-sJjuxs a solid 
crust significantly affects the size of the instability windowEllcSHa but has 
been suggested to have no effect for the evolution of young neutron stars. The 
size of the window is shown in Fig. 3, where we discuss the influence of the 
r-mode instability in accreting old neutron stars. 

Strange stars, if they exist, might consist of the most stable known form of 
matter, which is a mixture of strange, up and down quarks. There are specu- 
lations that a considerable number of the observed compact stars fall into this 
category. Both shear and bulk viscosity are significantly differertLfor strange 
stars, which has a significant influence on the instability windowtJ. Bulk vis- 
cosity is dominant for temperatures between 10^-10^" K and as a result the 
instability window is considerably smaller and moves to lower temperatures. 
As we will show later, this has a significant effect for Low-Mass X-Ray Bi- 
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naries (LMXB)li. In addition 



frequencies, see Fig. 7 inQ. 



a new instability window appears for higher 



In specific EOS it is expected that the central parts-oLthe neutron star 



the bulk viscosity due to non-leptonic weak interactions involving hyperons 
will be-considerably stronger. In a more detailed calculation Lindblom and 
OwenE^I have verified this result and have shown that the instability window is 
considerably smaller. They suggested that for young neutron stars, containing 
a considerable amount of hyperons in their cores, the r-mode instability will 
be suppressed before it significantly affects the angular momentum of the star. 

3.2 Evolution of nascent neutron stars 

Neutron stars are born after the gravitational collapse of the neutron-star 
progenitor (massive star or accreting white dwarf). Conservation of angu- 
lar momentum suggests that they should rotate almost at break-up speed, 
i.e. near their Kepler limit. Moreover, they should be quite hot with tem- 
peratures of the order of 10^^ K. These nascent neutron stars cool down by 
neutrino emission within a few seconds to temperatures below 10^" K with- 
out considerable loss of angular momentum. According to Fig. 2 the star 
then enters the instability window and the amplitude of the r-mode starts 
to grow exponentially. After a few minutes, nonlinear effects take over and 
saturate the highly unstable mode. When saturation occurs the star enters 
a nearly steady state, in which the amplitude of the mode remains essen- 
tially unchanged and the angular momentum of the star is radiated away 
while it gradually cools down. After a year or so the angular velocity and 
the temperature of the star are low enough that viscous damping can start to 
dominate over the r-mode instability. Consequently, the star moves out of the 
instability window and the phase of angular momentum loss to gravitational 
radiation ends. The star keeps cooling down slowly while its angular velocity 
remains practically constant. In the case of the best known pulsar, the Crab 
pulsar, this evolutionary scenario has provided some impressive results. The 
Crab pulsar has an observed period of 33 s, but 1000 years ago when it was 
born, it had an estimated period of 19 s, which is exactly the period that the 
evolutionary model predicts: Initially born with a period of 1-2 ms, the Crab 
pulsar has radiated away about 95% of its angular momentum in the form of 
gravitational waves within a year, thereby slowing down to a period of ~ 19 s. 

For young strange stars, it has been shown that the. r-modc driven spin- 
evolution is quite different from the neutron star caseE3. In a young strange 
star, the r-mode undergoes short cycles of instability during the first few 
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cores contain a significant 




months, followed by a quasi-adiabatic phase, where the r-mode remains at 
a small and roughly constant amplitude for thousands of years. Another 
important feature from the neutron_atar case is that the r-mode in a strange 
star never grows to large amplitudesES. These results suggest that the r-modes 
in young strange stars emit a persistent gravitational-wave signal that should 
be detectable, with large-scale interferometers, given an observation time of a 
few monthsEj. 



3.3 Gravitational radiation emission 

The aforementioned 95% loss of angular momentum due to the r-mode in- 
stability is transformed into gravitational radiation, which is emitted at a 
frequency range of 50-1000 Hz. This large amount of angular momentum loss 
in gravitational waves makes the rotational instabilities one of the primary 
sources for the detection of gravitational waves with laser interferometric de- 
tectors. 

The frequency of the waves produced by an f— j= 2, m — —2 r-mode is 
/gw — 217/37r and the expected strain amplitude i^ 

h{t) = 10-24 a Mi.4 i??o (^1^) , (22) 

where D is the distance from the source. The signal to noise ratio for the 
above signal is 

SV ^ r \Hf)?df 

where h{f) is the Fourier transform of the signal and Sh{f) is the one-sided 
power spectral density of the detector noise. If the saturation amplitude 
is about unity, newly formed neutron stars out to about lOMpc should be 
detectable by LIGO II with narrow-.banding. This suggests that the Virgo 
cluster is likely to be out of reacheil. For the first phase of the instability, 
when the signal grows up to a saturation point, there are estimates predicting 
an upper signal to noise ratio of S/N < 10 for the LIGO II configuratiorcS. 



3.4 Accreting neutron stars 

The r-mode instability revived the interest in the gravitational radiation pro- 
ducing mechanism-Suggested originally by Papaloizou and Pringleo and later 
also by Wagonercj. The mechanism applies to old neutron stars which are 
spun up by accretion until they reach the instability window. For these stars 
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it was assumed that the loss of angular momentum via gravitational wave 
should balance the accretion torque. 

This_spenario, initially suggestejd for /-modes, has been revived for the 
r-modesEJ'Ca. It was soon realizedEll that the viscous heating of the neutron 
star will not allow the star to remain in the equilibrium state suggested by 
Wagoner. Instead, the runway reheating will drive the star further into the 
instability regime (path A-B in Fig. 3) and spin it down (path B-C). This 
process will last for a few months after which the star will follow the path (C- 
D-A), i.e. it will cool down and spin up again due to accretion. Eventually, 
it will reach the instability window again and repeat this cycle. This quiet 
phase lasts for about 10^ years, which means that an accreting neutron star 
spends only a tiny fraction of its accretion lifetime in the phase A-B-C, during 
which it emits gravitational radiation. Still, Fig. 3 shows three remarkable 
resultal3. First the observed LMXB seem to follow the pattern suggested here, 
second this mechanism sets a lower limit for the rotation period of the recycled 
pulsars and third it provides an explanation for the period clustering of the 
millisecond pulsars. The point A in Fig. 3 corresponds to the smallest possible 
period that the pulsar can achieve via accretion. Remarkably, the theoretical 
value is 1.5 ms while the fastest known pulsar rotates with a period of 1.56 ms. 
Pulsars following this accretion scenario should have periods bounded by the 
two points A and D, which correspond to 1.5 and 6.5 ms and this is true for 
the majority of the known millisecond pulsars. 

Although Wagoner's original scenario is not realized for normal neutron 
stars, it applies to strange stars. The instability window is quite different here 
and it has been found that unstable r-modes affect strange stars in a way that 
is quite distinct from the neutron star case. For accreting strange stars, the 
onset of the r-mode instability does not lead to-Lhe thermo-gravitational run- 
away that is likely to occur in neutron stara23. Instead, the strange star 
evolves towards a quasi-equilibrium state on a timescale of about a year. The 
reason is that the instability window moves towards lower temperatures and 
the corresponding point A is on the left side of the window. When the star be- 
comes unstable, the runaway heating pushes the star to higher temperatures 
i.e. outside the window (ordinary neutron stars are pushed further inside). 
The star becomes stable again and soon cools down and "hits" the instability 
line again and the process continues for as long as the accretion lastsE3. This 
is shown in Fig. 4. The upper panel shows the time evolution of the mode 
amplitude a. After a few initial peaks, which are of the order of a ~ 10"**, the 
amplitude settles to a constant level of about 10~^. During this "stable" phase 
the mode heating remains also constant (middle panel) while the rotational 
period of the star increases because of the continuous loss of angular momen- 



florence: submitted to World Scientific on February 7, 2008 12 



0,0 1 ■ 1 ■ 1 ■ ' ■ ' ■ ' ■ ' ■ 1 

7.0 7.5 8.0 8.5 9.0 9.5 10.0 10.5 
Iog,„ (T/IK) 

Figure 3. The r-mode instability window relevant for old neutron stars. We show results for 
the simplest (crust-free) model (thin solid line), as well as for a star with a crust (thick solid 
line). We illustrate two typical r-mode cycles (for mode saturation amplitudes as =0.1 
and 1), resulting from thermo-gravitational runaway after the onset of instability. For 
comparison with observational data, we indicate the possible range of spin-periods inferred 
from current LMXB data (shaded box) as well as the observed periods and estimated upper 
limits of the temperature of some of the most rapidly spinning millisecond pulsars (short 
arrows) . 



turn in gravitational waves. This mechanism could explain the clustering of 
spin-frequencies inferred from kHz QPO data in LMXB. These results suggest 
that the r-modes in strange stars emit a persistent gravitational-wave signal 
that should be detectable with large-scale interferometers given an observa- 
tion time of a few months. If detected, these signals would provide unique 
evidence for the existence of strange stars, which would put useful constraints 
on the parameters of QCD. 

3.5 Magnetic Fields 

Neutron stars usually have strong magnetic fields, which should interact with 
a large amplitude pulsation mode. It has been suggested that a magnetic field 
should have a stroiig influence on the evolution of the r-mode instability and 
vice-versa. Spruitc3 suggested that the r-mode should generate some kind of 
differential rotation. In this case the magnetic field winds up until it reaches 
a critical point where it becomes unstable due to buoyancy forces. When this 
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Figure 4. The evolution of the r-mode amplitude (upper frame), the temperature (middle 
frame) and the spin period (bottom frame) corresponding to the first 10^" s following the 
initial onset of r-mode instability in an accreting strange star. 



happens the star emits a considerable amount of electromagnetic radiation 
and perhaps gamma rays. As a consequence the neutron star will acquire an 
extremely strong magnetic field aud-becpme a magnetar. 

Rezzolla, Lamb and ShapiroE2rEili£3 suggested that the magnetic field 
might actually prevent the r-mode instability at all. This is possible if the 
non-linear evolution of the r-modes winds up the magnetic field of the neu- 
tron star in such a way that it drains energy away from the mode and thus 
suppresses it entirely. For magnetic fields of the order oi B > 10^° G, the ro- 
tation rate must be above 0.35 Qk in order for the mode to survive, while for 
stronger magnetic fields of the order of i? > 10^^ G, the r-mode cannot grow 
at all. Ho and LaiO discussed the effect of magnetic braking on the r-mode 
growth and showed that it is important only for extremely strong magnetic 
fields B > 10^"* G. A more systematic study of the effect of the rnagnetic field 
on the r-modes was recently published by Morsink and Rezaniao. 

3.6 Nonlinear Calculations 

Many of the results that we presented till now were based on the assump- 
tion that the r-modes reach an amplitude of order unity and then saturate 
due to nonlinear processes. The way that the nonlinear effects influence the 
saturation of a mode is not known and there are questions which need to be 
answered. For example, does the coupling of the r-mode to other modes allow 
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the mode to grow to unit amplitude? Does the growth of the mode influence 
an initially uniformly rotating star in such a way that it becomes differentially 
rotating? We have seen that this is an important assumption for the influence 
of the magnetic field. 

Nonlinear studies of pulsation modes suffer from the immense computer 
power needed and both approaches that ate will discuss are limited by the 
computer resources. Stergioulas and FontE3 evolve a model of a rapidly ro- 
tating relativistic neutron star using the nonlinear hydrodynamical equations 
without evolving the spacetime (Cowling approximation). In this way one 
cannot see the actual growth of the mode, but by evolving a sequence of per- 
turbed configurations with increasing initial mode amplitude a, one can study 
at least the energy leak from the r-mode into other modes on a dynamical 
timescale. Though some other modes have been excited, no mode coupling 
has been observed at amplitudes of order unity at the resolution used in these 
simulations. 

These results were pjiZEiified by another nonlinear calculation by Lindblom, 
Tohline and VallisneriEjEj. They worked in Newtonian theory and used a 
post-Newtonian radiation-reaction force to drive the r-mode unstable. In 
order to shorten the growth time, they increased the driving radiation reaction 
force by a factor of 4500. In this way, they could observe the growth of the 
mode from an amplitude-p = 0.1 to a = 2. These results verified earlier 
approximate calculationsc^l of the loss of angular momentum of the star. They 
have not observed significant mode coupling either but at an amplitude a ~ 3.4 
they found nonlinear saturation, which they argue was due to the creation of 
shocks associated with the breaking of surface waves on the star. In both 
calculations it was observed that the r-xaodes generated, differential rotation. 

In two recent papers, Schenk at alE3 and MorsinkO developed a formal- 
ism to study the nonlinear interaction of modes in rotating Newtonian stars. 
Subsequently they usal this formalism to study the saturation amplitude of 
the r-mode instabilityEiJ. Their results suggest that the saturation amplitude 
is extremely small and never reaches amplitudes of order unity. Instead they 
showed that the mode saturates at an amplitude which is about four orders 
of magnitude smaller than that suggested by Lindblom et al.E3. This obvious 
disagreement raaidifi a consequence of the limited resolution in the nonlinear 
time evolutionsE^E3. |_. 

Recent estimations by Wagoneil£2l suggest, that even if one takes into 
account the assumption that the r-mode amplitude will never reach values 
of the order 1, the r-mode instability in accreting neutron stars will still be 
a very good source for gravitational jsiaves. This result is similar to the one 
presented earlier by Andersson et al.Ej for accreting strange stars. 
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4 Instabilities of rotating stars in GR 



The theory of non-radial perturbations of relativistic stars has been a field of 
intensive study for more than three-jdecades, beginning with the pioneering 
paper of Thorne and Campolattarolj. 

Tja^e authors focused on perturbations of non-rotating stars while 
Hartleys laid the foundations for computing rotating relativistic stellar models. 
He also devised a way of modelling slowly rotating stars. This slow-rotation 
approximation is still widely used since the problem becomes one dimensional 
and therefore much simpler than the two dimensional case of rapidly rotating 
and strongly deformed stars. 

Investigating the axisymmetric perturbations, Chandrasckhar and 
FerrariEj showed how rotation induces coupling of the polar and axial jnodes, 
which are decoupled in the non-rotating case. Soon after, KojimaEj pre- 
sented the first complete derivation of the coupled polar and axial perturba- 
tion equations. Rotation does not only couple the polar and axial modes, 
it also removes the {21 + l)-fold degeneracy of a mode with angular index I 
in a non-rotating star. To first order, this_Z£eman-like frequency splitting is 
linear in the rotation parameter. KojimaE3'Lj computed the splitting for the 
relativistic /-mode for various sequences of polytropic models. Yoshida and 
Kojimao compared the frequencies and their rotational corrections for / and 
p-modes in the relativistic Cowling approximation with the fully relativistic 
perturbation calculation. They concluded that the Cowling approximation is 
remarkably good and becomes better the more relativistic the stellar models 
are. 

The "run" on the r-modes started, when AnderssonEJ found that they are 
CFS-unstable for any rotation rate. In a first attempt to compute the growth 
time withitua fully relativistic framework, he used the axial equation derived 
by Kojimao, without the coupling to the polar equations. Although he found 
the right sign of the imaginary part of the complex r-mode frequencies, he was. 
not able to recover the 51^^^'+^) behavior (see eqn(19)). Ruoff and KokkotasEll 
pointed out that Kojima's equation is not adequate for computing r-modes, 
for it is derived by transforming the original axial field equations, which are 
first order in time (or in the eigenfrequency a), into a second order wave 
equation. In these manipulations, the eigenvalue a appears quadratically, 
and one obtains second order rotational corrections, which, however, were 
discarded. Since the r-mode frequency a is proportional to the rotation rate 
fl, it is clear that one has to keep the fl^ terms in the wave equation. 

Taking care of the proportionality between a and fl, KojimaO derived an 
approximate equation in the so-called low-frequency approximation, where 
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the gravitational radiation reaction is ignored. This equation is now an eigen- 
value equation, which is linear in a and contains only rotational correction 
terms linear in fi. Kojima immediately realized that the relativistic treat- 
ment introduced a new effect, not present in the Newtonian case, namely 
the appearance of a continuous spectrum, arising from the relativistic frame 
dragging, which introduces an effective r-dependcnt oscillation frequency. The 
presence of this continuous spectrum in Kojima's equation has been proven 
in a mathematically rigorous way by Beyer and Kokkotaalj. 

For quite a while, it was not clear whether this continuous spectrum 
was merely an artefact of the low-frequency approximation, which neglects 
the radiation reaction, or whether it occurred because of the neglect of the 
coupling to the polar equations, or a feature of the slow-rotation approxima- 
tion in general. Lockitch, Andersson and FriedmanO argued that Kojima's 
equation should only be valid for non-barotropic stars, and showed that for 
uniform density models, one can find regular mode solutions, representing the 
relativistic r-modes. _ _ 

Ruoff and KokkotasO and YoshidaO extended this study to polytropic 
models and found the somewhat surprising results that Kojima's equations 
yields r-mode solution only for a very restricted class of polytropic models, 
namely for those which small polytropic index n. For large n, no cegular 
solutions could be fouad.. Going one step further Ruoff and KokkotasE3 and 
Yoshida and FutamaseO included the radiation reaction, but still found that 
the basic properties of Kojima's equations do not change and its singular 
structure due to the continuous spectrum cannot be avoided. 

Kojima and Hosonumalla showed that if third-order rotational effects are 
added to the original Kojima's equation, it becomes a fourth order ordinary 
differential equation, which has no singular properties if the Schwarzschild 
discriminant associated with the buoyant force does not vanish inside the 
star. By solving a simpUfied toy version of the extended Kojima's equation, 
Lockitch and AnderssonLJ showed that the previously singular solutions can 
be regularized when the additional fourth-order term is taken into account. 

Yoshida and LeeEll and Ruoff, Stavridis and KokkotasES considered the 
fully coupled system in the Cowling approximation by including only first 
order corrections. The difference to the previous investigations is that the 
computation of the mode frequencies and eigenfunctions is not done in a 
perturbative method, in which Q. is regarded as a small expansion parameter. 
Instead the mode is a solution of the complete coupled system, with I ranging 
from ^ = m to infinity. For the numerical computation, this system has to 
be truncated at some cut-off value Zmax which has to be chosen such that the 
eigenfrequencies and eigenfunctions are well converged as Zmax is increased. 
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Yoshida and LeeEZI focused on non-barotropic models taken from McDer- 
mott, van Horn and Hansedlj, which have a finite r-dependent temperature. 
For all models under consideration, they found regular r-mode solutions, and 
all of them were lying in the regime of the continuous spectrum, which was 
forbidden from Xpjima's equation. 

Ruoff et al.ESl considered barotropic as well as non-barotropic polytropic 
stellar models and demonstrated that in both cases, one can always find r- 
modes. They discussed how the truncation of the coupled system of equations 
at some cut-off value /max affects the continuous spectrum, which for the cou- 
pled equations consists not only of a single band, but of several either discon- 
nected or overlapping patches. They found that some of the inertial modes, 
in particular the r-mode, can actually exist inside the continuous spectrum. 
However, for highly relativistic stellar models, the continuous spectrum can 
still make some of the inertial modes, which exist in less relativistic stars, 
disappear. 

In the remainder of this review, we will outline the various approximations 
and the results they are leading to. We will use geometrical units with G = 
c= 1. 



^.1 The slow-rotation formalism 



The metric describing a slowly rotating neutron star in spherical coordinates 
(i, r, e, (/)) is 



V 



2v 





e^^ 
r2 



'sin^6l 



-wr^ sin^ ( 



sin^ 6 



I 



(24) 



where v, A and the "frame dragging" uj are functions of the radial coordinate r 
only. With the neutron star matter described by a perfect fluid with pressure 
p, energy density e, and 4-velocity 

£-''(1,0,0,17), (25) 

the Einstein equations together with an equation of state p — p{e) yield the 
well known TOV equations plus an extra equation for the frame dragging, 
which to linear order is given by 



w" - Anre^^ip + e 



where 



4^ 

— I ZD' 



167re^^ (p + e) - 



0, 



(26) 
(27) 
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represents the angular velocity of the fluid relative to the local inertial frame. 

We assume the oscillations to be adiabatic, so that the relation between 
the Eulerian pressure perturbation dp and energy density perturbation Se is 
given by 

= (28) 

where Fi represents the adiabatic index of the perturbed configuration, T is 
the background adiabatic index 

r = ^^, (29) 

p e 

and f is the radial component of the fluid displacement vector If we write 
Fi as 

r.-^^ft). (30) 



p \de 



we can write relation ( P8|) as 



Sp= (^^) <5e-rA (31) 

with the Schwarzschild discriminant As given by equation (1). 

The metric perturbations h^i, are functions of all four spacetime variables 

h^, = h^,{t,r,e,^) (32) 

To eliminate the angular dependence, one expands the perturbations into 
tensorial harmonics, which are a generalization of the well-known spherical 
harmonics Yim = Yim.{0, (j))- One thus obtains a set of equations for the coef- 
flcients, which only depend on t and r and the parameters I and m. In the 
non-rotating case, these equations are independent of m and decoupled with 
respect to I. Furthermore, they fall into two independent sets, which show 
opposite behavior under parity transformation. The polar or even parity per- 
turbations transform as (—1)', whereas the axial or even parity perturbations 
as (—1)'+^. Rotation induces a coupling between the different sets, through 
flrst order rotational terms, which couple the set of polar/axial equations of 
order I to the set of axial/polar equations of order I ± 1. Furthermore, the 
equations now depend explicitly on m. 

In deriving the perturbation equations, one has to make use of the well- 
known gauge freedom. In the non-rotating case, the commonly used gauge is 
the Regge- Wheeler gauge. However, in the rotating case, it is more convenient 
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to use a diflFerent gauge, theJBCL gauge, which was first introduced by Bat- 
tiston, Cazzola and LucaronicJ. A convenient way to set up this gauge and to 
derive the perturbatdnn-equations for slowly rotating relativistic stars is to use 
the ADM-formalisirH'E2l. The ADM formalism provides a natural framework 
to derive the perturbation equations for rotating stars, as it clearly separates 
the equations into dynamical evolution equations and the set of time indepen- 
dent constraints. As a consequence the equations do not contain any mixed 
time and space derivatives and in the BCL gauge, they do not contain any 
second order derivatives at all. This makes them particular suitable for the 
numerical time evolution, and they can also be easily cast into an eigenvalue 
problem. 

In the ADM formalism, the gauge is chosen by giving prescriptions for the 
lapse function a and the shift vector As we are dealing with perturbations, 
a and are the perturbations of lapse and shift. In the BCL-gauge, one sets 
a = and is chosen such that the evolution equations for the angular 
perturbations /i^e, h^^ and h^^ become 

dthab = ^. a, 6 ={0,0}. (33) 

With hah = Oj the expansion of the metric perturbation is 

a = 0, (34) 

1.771 

^polar^Y^l ^ 1,™, (36) 



Lm V * 



for the polar part and 



^a.^al = ^ ( Q, -^3'™ ^ , V^"' siu Ode^ Yl^ , (37) 




-W"-?^ W™ sm0de 



^r"' = EU 0''"^ 1^"- (38) 


for the axial part. The asterisks stand for symmetric components. The pertur- 
bations of the extrinsic curvature are all non-zero and denoted by Ki through 
Kq. Of these, K3 and Kq represent the axial perturbations. Finally, the fluid 
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perturbations are decomposed according to 



5uf 



"Aal w \ ^ ( n „,lm ^4> „,lm ■ nfi \ \r 

— —e ^lU, — -U3 smtiao j Yim, 



l,7n 



l,7n 



1 



From Eq. (|28|), we obtain the relation 



P 



Tip 



(39) 

(40) 
(41) 

(42) 
(43) 



The requirement (Q) leads to the following three conditions for the shift 
components 



ylm ^ j^lrn^ 

Trim Ty^lm 

^2 — ^6 ^ 

'±1 



(44) 

(45) 
(46) 



where the operator is defined by its action on a perturbation variable 

plm 



<l — lm 



(; + 2)Q,+i„p'+i'", 



with 



elm 



' (/ - m){l + m) 
(2;-l)(2; + l)' 



(47) 



(48) 



These requirements completely specify the gauge, and we obtain three evolu- 
tion equations for the metric perturbations 6*3, V3 and V4, as well as for the 
six extrinsic curvature components Ki to K^,. From the perturbed energy- 
momentum conservation 5 {T^^^.^ = 0, we obtain the evolution equations for 
the enthalpy pertubations H and the velocity perturbations ui. In the non- 
barotropic case, we can derive the evolution equation for 1^ from the relation 
(p8|). In addition to the evolution equations, the Einstein equations yield four 
constraint equations which do not contain any time derivatives. The complete 
set of equations is rather lengthy and can be found inl£3. 
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4-. 2 Axial Equations 

We will first focus on the purely axial case. Neglecting the coupling to the 
polar equations, we have four evolution equations for the metric variable V4, 
the two axial components K:^ and Kq and the axial component of the 4- velocity 
U3 (in the following we omit the indices I and m) 



9 , A _ 1(1 + 1) -2^^ , 2im 



(49) 



- + imc^j Kf, = Vi - j^^^ [u:'K^ - 167r(r! - c^)(p + 6)^3] , (51) 

d \ 2mi{Q. -uj) ^ 

' imQ U3 ^ — , , — [U3 - Ke) , (52) 



^dt J " 1(1 + 1) 

together with the momentum constraint equation 

16Hp + e)u, = K^ + - fcll^i^e - |^e--V^, . (53) 
r /(/ + 1) 

In the literature, the two axial metric components are normally denoted by 
/lo and hi. In the above equations, V4 is proportional to hi and Kq is propor- 
tional to /iQ. In the non-rotating case M3 = 0, and the remaining equations 
( ^ ) - (^ are equivalent to a single wave equation for V4, which in the ex- 
terior reduces to the well-known Regge- Wheeler equation. For neutron stars, 
these equations describe thc^axial w-modes, which have been first studied by 
Chandrasekhar and FerrariEj for ultra-relativistic models, where they become 
long lived "trapped modes" , and by KokkotasEJ for less compact models. 

In the rotating case, fluid motion is possible, and a new family of modes 
is expected to appear, the unstable r-modes. For ultra-compact and rapidly 
rotating stars, an ergoregion can appear, which is able to drive some of the 
w-modes uiiS|table. This was known only for scalar and electromagnetic pefc 
turbations and has been investigated in detail by Commins and Schuta53 
for scalar fields on a stellar background. In the following we will discuss it for 
the gravitational perturbations themselves. 

4. 3 The w-mode instability 

For modes mainly associated with fluid motion, it is the rotation of the star 
itself which is responsible for dragging a counter-rotating mode forward. For 
w-modes, which are predominantly spacetime oscillations it is not the stellar 



florence: submitted to World Scientific on February 7, 2008 22 



rotation itself, but the associated frame dragging, that influences the behavior 
of the w-modes. In the non-rotating case, the axial w-modes do not couple 
to any fluid oscillations at all. If the star is set into rotation, this remains 



only approximately true as the fluid motion is induced by Eq. (52), and in 
principle also through the coupling to the polar equations. 

But if one still regards the w-mode oscillations as being essentially in- 
dependent of the fluid, then it becomes clear that the only way to make a 
u)-mode unstable is to create an ergosphere, i.e. a region where no timelike 
counter-rotating geodesies exist any more, or in other words, a region where 
the dragging is so strong that any timelike backwards moving trajectory gets 
dragged forward. This is associated with a change in the,aign of gtt- 

To test the existence of an ergosphere we write gu astj 

gtt = -e^" + u?r^ sin^ Q (54) 

and restrict ourselves to the equatorial plane B — -k 12. An ergosphere starts 
to appear whenever gtt > 0, or equivalently when 

<ruj . (55) 

Models allowing ergospheres have been constructed in the_case of rapidly ro- 
tating uniform density models by Butterworth and IpscrEll and for toroidal 
polytropic configurations by Komatsu, Eriguchi and HachisuH. It is clear 
that for any realistic equation of state, it is very unlikely to construct stable 
stellar models which would exhibit an ergosphere, nevertheless recent results 
suggest the existense of ultra compact quark starsO Therefore, according 
to the present observational data and the known types of EOS, one would 
not expect the ergosphere instability of the w-modes to be of any significant 
astrophysical relevance. The compactness needed for a star to develop ergo- 
sphere is extremely high i.e. R/M ~ 2.26-2.4. An interesting feature of this 
instability is that it can set in even for quite low rotation rates. For example 
for a star with R/M = 2.26 the first w-mode becomes unstable at a rotation 
rate of 0.19 ^Ik- For R/M = 2.4 this occurs for 0.83 17 Given that a star 
is becoming unstable due to the existence of an ergosphere then there is no 
direct dissipation due to shear or bulk viscosity, since the w-modes practically 
do not couple to the fiuid motions! This is a intriguing result since it suggests 
that the instability will grow and the star will radiate away in gravitational 
waves all its available angular momentum before it becomes stable againcl. 

4-. 4 The r-mode instability 

In contrast to the lu-modes, which are predominantly spacetime oscillations, 
the r-modes are characterized by strong axial fiuid currents. Even though, it is 
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Figure 5. Appearance of the ergoregion as a function of the rotation rate e = Q/Qj^. For 
the stellar model with R/M = 2.4, it first appears at r = 0.2R for £ = 0.15 and expands as 
e is increased. The vertical lines represent the rotation rates, at which the first w-mode of 
the respective models becomes unstable. 



only through the metric perturbations, that the r-modes are driven unstable, 
it is instructive to study their features in various approximations, where some 
or all metric perturbations are neglected. 

The most radical way is to neglect all metric perturbations. For modes 
with small damping, such as the polar p-modes for example, this so-called 
relativistic Cowling approximation yields quite accurate results. In Newtonian 
theory, it also yields very reliable results for the r-modes. Let us therefore 
apply it to the relativistic case as well. Neglecting all metric perturbation in 
our set of axial evolution equations (|4^) - (js^), the only remaining equation 
is that for the fluid variable M3, which can be written as 

1-3 = -^m (n - U3 . (56) 

From this equation wc can immediately deduce that the fluid layers are decou- 
pled from each other, which means that each layer has its own real oscillation 
frequency given by 

a = -m il - — — 



(57) 
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R/M = 12% 
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Figure 6. The real part of the frequency of the three ui-mode together with the r-mode of 
the stellar model with R/M = 2.28. The toi-mode crosses the zero line at e = 0.32. 



The whole frequency range is determined by the values of w at the center and 
at the stellar surface, and is the larger the more relativistic the stellar model 
is. In the Newtonian limit (tu ^ O), Eq. ( ^7| ) yields the well known r-mode 
frequency 



^^ = -mn[l-j^^y (58) 

In the relativistic case, the presence of the frame dragging uj destroys the 
occurrence of a single mode frequency and gives rise to a continuous spectrum, 
at least to this order of approximation. 

This somewhat peculiar result might suggest that the Cowling approxi- 
mation fails in the relativistic case, and we should rather keep some of the 
metric perturbations. To assess which ones to keep, we consider the non- 
rotating case. There ii^ = 0, i.e. we have no oscillatory motion. We can have, 
however, have a stationary fluid current with a nonzero associated metric per- 
turbation, which we can easily find by setting in our evolution equations ( ^9| ) 
- ( ^2|) all the time derivatives to zero (together with the rotational terms). 
From equations (^o|) and (|5l|), we immediately see that V4 has to vanish and 
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from Eq. (49) we obtain 



X' 



„2A 



^3 = . 



(59) 



Moreover, we have the constraint (|53|), which we can combine with the above 
equation to yield a single equation for Kq. To this end, we differentiate ( p9| ) 
with respect to r and use the constraint ( |53| ) to eliminate K^. This gives us an 
elliptic second order equation for Kq with a source term containing the fluid 
velocity U3 . Instead of using Kq , it is more convenient to write this equation 
in terms of the variable hn ~ e'^^^Kp. 



e~^^h'^-A-Kr{p + e)h'^ + 



^ AM l{l + 1) 
87r(p + e) H j ^ — 



/lo = l&TTe''-^{p + e)u3,. 

(60) 

Although this equation is derived under the assumption of stationary per- 
turbations, we can also use it for slowly rotating starts and low-frequency 
oscillations. In this case we no longer assume — 0, but we use the evolu- 
tion equation ( |5^ ) to govern the time dependence of M3. 

It is obvious that the above equation (|60| ) together with the evolution 
equation (|2|) is qualitatively different from the Cowling approximation (|56|). 
In the latter case, the fluid layers are completely decoupled and each one 
oscillates with its own frequency. With the inclusion of /iq, we have an elliptic 
Poisson-like equation, with the fluid acting as a source. As /iq couples back 
into the evolution equation (|5^) , it connects the fluid cells and might give rise 
to some coherent motion. 

To find out whether we can find mode solutions, we assume a harmonic 
time dependence of U3, i.e. M3 exp(icri). With this ansatz, we can combine 
into a single equation, which can be written as 



equations (52) and (F 



[a — w) 



-t- 167r(j3 + e)ahQ — , 
where we have introduced the normalized frequency 



1)1 



—) 



and 



The relation between U3 and ho becomes 

w 

U3 = e 



A- 



(61) 



(62) 



(63) 



(64) 
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Figure 7. The boundaries of the continuous spectrum vja and ros together with the r-mode 
frequency a as a function of the polytropic index n for stellar models with compactness of 
M/ R = 0.2. The case n = corresponds to a uniform density model. 



Eq. (|6l|), which is Kojima's master equationEl, is a singular eigenvalue prob- 
lem, because the factor a ~ w becomes zero if a falls into the range 

Wc < oi < Wg ■ (65) 

In this case, the fluid perturbation diverges urdess vanishes at the singular 
point, which one can prove not to be the caseES. If we translate the singular 
range for a back into the corresponding range for the eigenfrequency ct, we 
see that it is exactly the same as the continuous frequency spectrum we found 



in the Cowling approximation given by Eq. (57) 



One can show that mode solutions can only exist for a < lE3. If a lies 
in the range < ct < I, the eigenvalue equation is no longer singular, as 
the singular point now lies outside the star, where U3 is zero. Only for such 
values of a, regular mode solutions can be found. _ 
For uniform density models one can indeed find regular mode solutionsO. 
However for polytropic equatiopSpXif state, this is only possible for a very 
restricted range of stellar modelsEjO. The existence of regular mode solutions 
strongly depends on the polytropic index n and the compactness of the stellar 
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Period in ms 



Figure 8. Double logarithmic plot of the growth time for uniform density models with 
compactness M/R = 0.1, 0.2. Also included is the growth time deduced from formula (po|), 
which agree well with the numerical results. 



model under consideration. The general picture is that the smaller n is, the 
larger is the compactness range where one can find physical mode solutions. 
For a given polytropic index n, one usually finds physical mode solutions for 
models with a small M/R ratio. As the compactness is increased, i.e. as the 
models become more relativistic, the mode frequency a decreases and starts 
approaching Wg ■ Eventually it crosses this point and migrates inside the range 
of the continuous spectrum, thus becoming unphysical, and no r-mode exists 
any more. This is shown in Fig. ^ 

It has been argued that the singular structure which is associated with 
the continuous spectrum is an artefact of neglecting the radiation reaction. 
Both in the Cowling and low-frequency approximations, the eigenfrequency 
a is always real, whereas it actually should be complex since the emission of 
gravitational waves should either damp the modes or make them grow. If the 
frequency is complex, the singular point is removed from the real axis, and 
the differential equation is no longer singular. 

That the inclusion of the radiatinn reaction does not resolve this issue has 
been shown by Ruoff and KokkotasEZl by solving the full set of axial equations 
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and by Yoshida and Futamase0 by imposing a near-zone boundary condition 
for Kojima's equations instead of the asymptotic flatness condition. 

The conversion of the evolution system - (^) into a time independent 
form is straightforward. Assuming a harmonic time dependence e*' 
deduce the following two ordinary differential equations 



Vi 



1(1 + 1) 



[lu' — 1Qt:vj{p + e)u^) + i (cr + muj) Kq 



„2A 



X' 



together with the two algebraic relations 
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1(1 + 1) 



2imuj' 



(a + moj) Vi 



2mw 



2mw - 



1(1 



1) 



we can 

(66) 
(67) 

(68) 
(69) 



1(1 + 1)((7 + mn)' 

Again, Eq. (|6^) becomes singular if a is real and falls in the range of the 
continuous spectrum. For real a the denominator a + muj in Eq. ( |68|) can also 
vanish and the equation becomes ill defined. However, as a is expected to be 
a complex frequency, both equations for K3 and 7/3 should remain regular. 

For uniform density models, the low-frequency approximation always 
yielded r-mode solutions. The inclusion of the radiation reaction should then 
yield the associated growth times. For stellar models, which are not too com- 
pact, the fully relativistic calculation can indeed reproduce the Newtonian 
value given for Z = m = 2 byE3 



TpN = 22 



1AM, 



M 



10 km 

R 



P 
1 ms 



(70) 



with a remarkably good agreement. For very compact models, however, the 
growth times found in the fully relativistic case are in general larger than 
suggested by formula (^0|). In Fig. |^ we show the growth times for two uniform 
density models with M/R = 0.1 and 0.2 as function of the rotation period. 

When switching to polytropes, we still find the same pictures as in the 
low-frequency approximation. As long as the mode is well outside the contin- 
uous spectrum, we find a non- vanishing imaginary part, but as soon as the 
mode approaches the continuous spectrum, the imaginary part starts to de- 
crease and drops to zero at the point where the real part enters the continuous 
spectrum. This is depicted in Fig. |9| where the dependence of the complex 
eigenfrequency cr on the polytropic index n is shown. For small n, the real 
part of a lies well outside the continuous spectrum, but as n increases, the 



florence: submitted to World Scientific on February 7, 2008 29 



0.036 



1.5x10" 



-8 




Re(a) 



0.029 - 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



,0.0 



n 



Figure 9. Behavior of the real and imaginary part of the r-mode frequency cr as a function of 
the polytropic index n for models with M/R = 0.2 and rotation period of 1.0ms. The scale 

on the left is for the real part, the scale on the right for the imaginary part of cr. For small 
n, the real part Rc((t) lies outside the range of the continuous spectrum. For increasing n, 
Rc((t) approaches the lower boundary and the imaginary part lm(a) eventually starts to 
decrease. For n = 0.843, 5R((t) reaches the lower boundary of the continuous spectrum and 
Im((T) goes to zero. For larger n, no mode can be found any more. 



continuous spectrum expands and approaches the real part of a. The imag- 
inary part first grows with n, but as the continuous spectrum expands and 
approaches the real part of cr, it bends down and goes to zero very rapidly. 
For n = 0.843, is becomes zero and for larger n, a regular mode no longer 
exists. 

4-5 Inclusion of the coupling 

The above investigations have shown that purely axial modes do not necessar- 
ily exist in relativistic non-barotropic stars. So far we have assumed that the 
coupling to the polar equations can be neglected. However, as the results of 
the previous sections have shown, this might lead to very unexpected results. 
To assess, whether the existence of the r-modes indeed depends on the stellar 
model, one has to go one step further and include the coupling to the polar 
equations. As the complete set of equations is quite complicated, it is useful 
to study the effect of the coupling within some further approximations. 

Although the Cowling approximation did not give any useful results in 
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the purely axial case, it nevertheless yields the key to solving the puzzling 
results of the previous sections. 

The first order evolution equations in the Cowling approximation are 
derived from the perturbed energy-romncntum conservation law S (T'"'.^) = 0, 



and read in the non-barotropic cas' 



21/-2A 



{dt + imn) H = e 



{dt + imn) ui = H' + 



C 



{dt + imn) U2 — H 



V_ 

riP 



;eB 

2v' - y 

2imvje^^'^''H 



2 

r 

v'ui 



Ui 



2x1(1 + ^) 



U2 



1(1 + 1) 



(imu2 + C^^u^ 



{dt + imn) U3 



2vj 



l{l + l) 



{imus — Cf^U2) 



(71) 

- B {imu2 + Cf^us) , (72) 
Bu,, (73) 
(74) 



y-.2 g 2A 

'W+T) 



j,2 g 2A 



B£f\i, 



{dt + imn) C = i/' 



1 e 



21/-2A 



Ul, 



with 



B = w' + 2tu (^i^' - i j , 
and the coupling operators 

^±lplrn ^ _ 1)Q,„P'-1™ -{1 + 2)Qi+i„,F'+l'", 

-{l + l)QimP'-^"' + lQi+i^P'+^"\ 

l{l + 2)Qi+i„,P 



-±1 ul 



Cf'P'"' ^ {l~l){l + l)Qi„,P- 



i-ln 



l + lm 



(75) 

(76) 

(77) 
(78) 
(79) 



As these operators couple equations of order I to equations of order ^ ± 1, the 
above set of equations represents a coupled system with / ranging from I = m 
to ^ = 00. 



4-.6 The Continuous Spectrum 

By assuming a harmonic time dependence, one can easily transform the above 
time dependent equations into an eigenvalue problem for the frequency a. 
Defining 
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the eigenvalue equations can write be written as a set of ODEs 

(s^ = As^ + bLu'-^ + B^u^ + B[u^+^ l^m,...,oo (81) 

together with the algebraic relations 

UW'^ + Uy/ + U[u^+'^ = SLs^-^ + SKs^ + S-'^s'+i, l = m,.. . ,oo,(82) 

where the letters A, B, U and S with their various indices represent 2x2 
matrices. The algebraic relations can be viewed as a matrix equation for two 
infinitely dimensional matrices U and S acting on the vectors 

u = (u'",u™+i,u™+^...,u',...)^, (83) 

s= (s™,s™+\s"+^...,s',...)^, (84) 

whose respective elements are the 2-vectors u' and sK Both U and S are 
tridiagonal block matrices, with each block given by the above 2x2 matrices. 
Now we can write Eq. (|8^ ) as 

^U"'w'' =^S"'s'' , / = m...oo, (85) 

I' I' 

which can be solved for by multiplying both sides with U""'^ 

= J2{U-^f' S'''"/' , l = m...oo. (86) 

I'l" 

Note that the matrix U is r-dependent since its elements C/j_ contain the 
function zu. For U to be invertible its determinant must not vanish. For 
certain values of ct, however, this happens. In the Newtonian limit m fl, 
the matrix U becomes independent of r, and one can easily show that U 
becomes singular only for a discrete set of frequencies a. These frequencies 
represent the solutions of the homogeneous part Uu = of Eq. (^5|). In the 
relativistic case, however, the zeroes of det U depend on the position inside 
the star, which means that each single Newtonian frequency will be spread 
out into a continuous band of frequencies, determined by the values of vu at 
the center and at the surface of the star, which we will denote by zuq and -cufi, 
respectively. In order to numerically perform the inversion and solve equation 
(^), we have to truncate the system at some value Ima-x ^ fn- To explicitly 
compute the ranges of the resulting continuous spectrum we need the matrices 
U'_, [/^ and which are given by 

f/i=-2ti7^g„„(^Jj), (87) 
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U[ = -2^^J±^Ql+^,r^ (liy (89) 

with 

^ 2r7iti7 . ^. 

and 

K = a + mn. (91) 

In the following examples, we choose m — 2. Let us start with the simplest 
case which is choosing Z„iax = 2. Then U — and det U = just yields 

S = 0, (92) 

which gives exactly the same continuous spectrum as following from Eq. (|57|). 

For Iniax ^ 3 

and det U = yields the two solutions 

-=(^2±^j- (94) 

Hence, the two ranges of the continuous spectrum are given by 

1.2319 Tuo < K < 1.2319 OTfl, (95) 

for the plus sign and 

- 0.2319 u7o<K< -0.2319 wr, (96) 
for the minus sign. For /„iax = 4, we obtain 



K 



1.4964 

0.4669 . (97) 
-0.7633 



Each increment of Zmax leads to an additional range of the continuous spec- 
trum. Depending on the respective values of vjq and tu^, the various ranges 
might as well overlap. This actually happens for the more relativistic stars 
since the variation in vj is much greater. 
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Figure 10. The continuous spectrum and the r-mode for m = 2 as functions of M/R for 
/max = 2 and imax = 3. For /max = 2, only the continuous spectrum exists. The relativistic 
r-mode appears for /max = 3 well outside the continuous spectrum as the latter has shifted 
away from its previous location for /max = 2. The dotted line represent the Newtonian 
r-mode value from Eq. (pSl). 



In the Newtonian limit tn — > f2, the continuous spectrum shrinks down 
to discrete vahies. Interestingly, these values n/ft, following from det U — 0, 
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Figure 11. Same as in Fig. 10 for m = 2. In addition we have plotted the Newtonian value 
for the r-mode according to Formula (js^ . The relativistic r-mode appears for (max = 3 well 
outside the continuous spectrum as the latter has shifted away from its previous location 
for I max — 2. Only the r-mode and the 22 and 22 modes can exist for the whole compactness 
range. 



are exactly the eigenvalues kq of the (generalized) r-modes of the Maclau- 
rin spheroids in the limit of low angular velocities, which were computed by 
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Lindblom and IpserEll and presented in their Table 1 . Our above values have 
been obtained for m = 2, but we can reproduce all their values for m = 
to m = 3. Note, however, that their I is related to our ^max by I = ^max — 1- 
Also note that these values are completely independent of any stellar structure 
parameters. 

In Figs. 10 and |ll| we show both the continuous spectrum and the eigen- 



modes as functions of the compactness M/R for different values of Zmax- The 
modes are for bamtropic stellar models, based on a polytropic equation of 
state with n — ILj. Starting for /max — 2 (upper panel of Fig. |l^), we find 
no mode, but only a band of frequencies given by Eq. (|57|). We have also 
included the Newtonian value for the r-mode as given by Eq. (^8|). One can 
clearly see that in the limit M/R the continuous spectrum shrinks to 
this value. If we include the coupling to the polar equations with I — 3, 
the situation changes drastically. As expected we now obtain two continuous 
frequency bands and a mode which lies below the Newtonian r-mode value, 
the further below the more compact the star is. Translating to positive fre- 
quencies, this means that this mode always has a higher frequency than the 
Newtonian r mode. Inspection of the eigenfunction reveals that the axial part 
is well described by a r''~^^ power law, i.e. 

which is exactly what one expects for the r-mode. Inclusion of higher I only 
gives minor frequency corrections, which tells us that this is indeed a physical 
inertial mode, the relativistic version of an to = 2 r-mode. According to our 
nomenclature for inertial modes, the r-mode is actually an if inertial mode. 
In Fig. In], we have also included the coupling up to Zmax = 4 and 5. Here, 
each inclusion of one additional I leads to another patch in the continuous 
spectrum, and to the appearance of some new inertial modes. Notice that 
the width of the uppermost and lowermost patches is very broad, whereas the 
middle patches around <j/^k ~ — 1 remain very narrow, even for large M/R. 
This is because the width of a patch is given by k {tur ~ '^0)1 which is therefore 
the smaller the closer k is to zero. For any given Zmaxj we find Imax ^ ^ new 
inertial modes. We label them with I = /max ^^nd n = 1, — to. In the 
Newtonian limit, they appear very close to the values found from dct U = 
for /max — 1- This is what to be expected from the Newtonian theory of 
inertial modes, since the eigenfrequency can be computed at lower order than 
the eigenfunction. This means that for the inertial modes appearing at 
/max = /, the associated eigenfrequencies can already be inferred from the 
solutions of det U = for /max — 1- 

From Fig. |o|, it is clear that the r-mode exists for the whole compactness 
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range, as it is always located outside the continuous spectrum. For ^max = 4 
(upper panel of Fig. ^l]), it is still outside the continuous spectrum, but for 
^max = 5 (lower panel of Fig. ^l|), it lies inside the continuous spectrum 
for M/R > 0.14. In this case, the eigenvalue equations cannot be solved 
any more, as U becomes singular. It is only through the evolution of the 
time dependent equations that one can see that the r-mode exists inside the 
continuous spectrum. This is not only the case for the r-mode, but also for 
other inertial modes. However, not all the modes are able to survive inside 
the continuous spectrum. This can be seen for the inertial modes if and in 
Fig. ^ For ^inax = 4, they both cease to exist when they reach the continuous 
spectrum. For Zmax = 5 the continuous spectrum has shifted such that the 
?2 mode remains always outside whereas the if mode still cannot penetrate 
the continuous spectrum. For larger Zmax, one can find the i| mode inside the 
continuous spectrum, but not the if mode. 

These results can be summarized as follows. As a result of the relativistic 
frame dragging, there always exists a continuous spectrum, which, however, 
depends very strongly on the number of coupled equations. For ^^ax = 2, the 
continuous spectrum is confined to a single connected region, with increasing 
width for more relativistic stellar models. If l^ax is increased the continuous 
spectrum splits into an increasing number of patches, which are disconnected 
for weakly relativistic stars but tend to overlap for more compact stellar mod- 
els. In the limit Zmax oo the continuous spectrum should completely cover 
the frequency range of the inertial modes. Nevertheless, there are inertial 
modes, in particular the r-mode, which can exist inside the continuous spec- 
trum. For very relativistic stars, however, the continuous spectrum is able to 
destroy some of the inertial modes. 

Caution should be used in order to draw any conclusions about the 
(non)existence of these modes for a rapidly rotating neutron star. The above 
results showed that the continuous spectrum is very sensitive to the number 
of equations that are coupled. So far, only first order rotational corrections 
have been taken into account, i.e. only the coupling from I to Inclusion 
of higher order corrections leads to coupling to higher I, which can modify 
the continuous spectrum in a significant way. It is therefore possible that 
higher order corrections affect the continuous spectrum in such a way that 
modes which do not exist in the first order analysis can now exist because the 
responsible patch could have moved elsewhere. 

Of course, another important step consists in including the metric pertur- 
bations. Lockitch et al.O included some (non-radiative) metric perturbations, 
but restricted their studies to a post Newtonian treatment. Still, the continu- 
ous spectrum will remain unaffected when metric variables are included, for it 
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is only the fluid equations which are responsible for its existence. However, by 
including the metric perturbations, the modes might be affected in such a way 
that they get pushed out of the continuous spectrum and we might find some 
of the modes which we cannot find in the Cowling approximation. That this 
might true can be seen from the purely axial I = m case, which we have dis- 
cussed in the previous sections. There, the Cowling approximation leads to a 
purely continuous spectrum without any mode solutions. In the low-frequency 
approximation, where one retains a certain metric component, one is lead to 
Kojima's master equation, which in certain cases admits r-mode solutions, 
whose frequencies lie outside the continuous spectrum. But the range of the 
continuous spectrum is still the same as in the Cowling approximation and 
even when all the metric perturbations were included, the qualitative picture 
remained the same. 

Nevertheless it is clear that the relativistic case is quite different from the 
Newtonian one, and further studies are necessary to gain a full understanding 
of the inertial modes of (rapidly) rotating relativistic stars. 
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